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Abstract 

Recent progresses in understanding the behavior of dense granular flows are presented. After presenting 
a bulk rheology of granular materials, I focus on the new developments to account for non-local effects, 
and on ongoing research concerning the surface rheology and the evolution of mechanical properties for 
heterogeneous systems. 

Resume 

Proprietees rheologiques des ecoulements granulaires denses. Les recents pro-grfe de comprehension 
des ecoulements granulaires denses sont passes en revue en se focalisant sur la rheologie classique des 
materiaux granulaires tous en rafhnant notamment sur les effets non-locaux, I’influence des conditions limites 
des ecoulements et revolution des proprietes d’ecoulements de systemes heterogenes. 
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1. Introduction 

Granular materials have been extensively studied in particular due to their large involvement in geo¬ 
physical events and in industrial processes. From powders to boulders avalanches, the granular systems are 
composed of individual grains interacting through solid contacts and they share some features. Compared 
to classical Brownian particles or molecules, their relatively large size (d > 100 /im) prevents thermal energy 
from significantly sustaining agitation or leading to an equilibrium state. Thus external mechanical energy 
must be supplied to shear or to deform the system. 

Depending on the amount of energy, three qualitatively different states can be distinguished [1]. When 
the grains have very little kinetic energy, the system is dense, the structure is dominated by force chains and 
the behavior is well described by the soil mechanics. On the other hand, if a lot of energy is brought to the 
grains, the system is dilute, grains interact only by binary collisions and the flow properties share common 
features with an agitated gas. 

However, due to the dissipative nature of these interactions, the loss of energy may decrease the kinetic 
energy leading rapidly to a denser system as often encountered in common granular flows. In such flows, 
grains have enduring long-lasting contacts with their neighbors and mostly interact through frictional and 
collisional forces. Despite the large interest for granular flows and the recent advances in modeling, a general 
theoretical framework for a granular rheology covering all the situations is far from been achieved. In this 
article, I will mainly focus on the rheological properties of spherical hard grains interacting through dry 
frictional contacts. One example is the flow on inclined plane as illustrated in figure Q where I split the 
problem of describing the flow into three numbered issues. The first one concerns the lowest-order laws 
able to model the velocity profile in the bulk. The second one deals with spatial correlations and their 


Email address: pierre.jopOsaint-gobain.com (Pierre Jop) 


Preprint submitted to C. R. Physique 


March 19, 2015 






Figure 1: Schematic of a flow on inclined plane illustrating three issues to describe the dense granular flows. First, what can 
predict a frictional rheology for the bulk flow? Second, when are the coupling across space important? Third, which interaction 
law governs the flow properties at boundaries? 


consequences on the rheology. The third one aims to describe the interactions between the flow and the 
boundaries. These issues will guide the discussion along this article. 

In the following, I will first present the bulk behavior of the granular matter and a rheological law (section 
[^. Since deviations from its predictions exist, I will then show different ways to include ”non-local” effects 
to extend the previous theoretical description (section]^. In section]^ I will make an overview of the surface 
properties issues that arise in real systems and the different models trying to incorporate them. In the last 
part (section , I will discuss the evolution of this model when dealing with other geometries of grains or 
interactions. 


2. Bulk rheology 

In this section, I will assume that the local friction coefficient does not depend on higher order derivatives 
of the velocity field than the local stain rate. For sake of simplicity, I will call this rheology a local rheology. 
On the contrary, if higher-order derivatives in space are involved or if other helds modifies the lowest-order 
behavior leading to a stronger coupling over distance, we cannot define a unique flow curve linking the shear 
stress to the strain rate (and to the pressure). I will call it a non-local rheology in this article and I will 
discuss this point in the section]^ 

2.1. Properties of steady dense flows 

I briefly describe here the basic features of steady dense granular flows in simple geometries. To have a 
complete overview of flow properties both in experiments and in numerical simulations, readers may refer 
to some recent reviews mw- 

From a rheological point of view, the simplest conhguration is the plane shear geometry. The granular 
layer is only confined by an external constant pressure P and sheared by driving one boundary either at 
constant velocity or at constant shear stress r. In steady state, simulations have shown that the velocity 
profile is almost linear except close to the boundaries [SJE]. 

Experimentally it is more convenient to study surface flows. The inclined plane configuration has the 
advantage to produce uniform flows, for suitable parameters, with only two control parameters: the inclina¬ 
tion of the plane (d) and the thickness of the flowing granular layer h. Above a thickness threshold (hstop(^)) 
the material flows and exhibits steady uniform flows in a narrow range of inclinations of few degrees (5° or 
10°) [TIISIIS]- For high enough thicknesses, the mean velocity scales with the thickness to the power 3/2 and 
is inversely proportional to the thickness of the deposit hstop(^)- This scaling is known to fail close to the 
threshold, for thin layers [IHIIII], I will discuss this point in section]^ These observations were identical in 
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numerical simulations [5] .When the grains are not flowing on a plane but on top of a pile, the convexity of 
the velocity profile changes presenting a fast flowing layer and a creeping flow deep in the pile |12) . In this 
configuration, the slope is selected by the system and the flow rate [13]. Finally, in confined conhgurations 
such as the Couette cylinder or the silo where the stress spatial distribution is highly heterogeneous, the 
flow exhibits shear localization close to their boundaries [5|. 

Some of these behaviors can be gathered in the same theoretical framework. In the following section, I 
will introduce the main ideas of this rheological model. 

2.2. Dimensional analysis 

Assuming a collection of rigid frictional spheres in a plane shear geometry, one can neglect the elasticity. 
Then using dimensional arguments, one can define only two dimensionless quantities: the ratio of shear 
stress and normal stress defining a friction coefficient fj, and an inertial parameter I Hi], also known as 
Savage number m or Coulomb number m- 


I = 

y/P/Pp' 


( 1 ) 


where 7 is the shear rate, d is the grain diameter, P is the constant confining pressure and Pp the grain density. 
This number can be viewed as the ratio of two time scales: The first one associated to the macroscopic shear 
rate is the time needed for a grain to overtake its neighbor (I/ 7 ). The second is associated to the microscopic 
rearrangement, it is the time required for a grain to fall behind its neighbor under the confining pressure 
{\Jd?pp/P). Decreasing /, the system gets closer to the quasi-static regime, whereas the kinetic regime is 
reached for high value of I. Having only two dimensionless parameters, the friction coefficient p.{I) should 
be only a function of I. This was indeed proved by discrete numerical simulations of plane shear [SlIS]. For 
the dense flow regime, corresponding to moderate inertial number (from O.OI to 0.5), the friction coefficient 
/i(/) and the volume fraction can be modeled accurately by the following phenomenological equations: 


M = Ms + and (j) = (j)ma.x - ( 2 ) 

io/t + i 

where the constants are material dependent, Iq ~ 0.3, « 0.2. Figure ([^left) shows the typical behavior 

of the friction coefficient p, and volume fraction (j) with the inertial number I. The friction coefficient starts 
then from a minimum value ps which coincides with the internal Coulomb friction of the material (typically 
Ps ~ tan21° for glass spheres), and saturates toward a value pg {pg « tan32°). One can observed that 
experimentally, steady homogeneous flows with higher values than I = 0.5 are difficult to achieved due 
to instabilities or transition to the kinetic regime (see section i 0111, thus the validity of Eq. at high 
I is questionable. The asymptotic limit pg is nevertheless required to account for the development of a 
gaseous phase in front of an avalanche [T^ or on top of a flow over a steep slope (tan0 > pg) [T7|[T3]. For 
practical applications, since the volume fraction remains close to its maximum value, one can often neglect 
its variations as a first order approximation and use only the evolution of p with I. 

It is worth noting that these laws were discovered in the plane shear simulations when the volume fraction 
was no longer set. Instead, the pressure is controlled and the volume fraction is free to adapt in the system 
with the evolution of other parameters. On the contrary if the volume fraction is fixed and the normal 
pressure measured, this will fixed the value of /. Then Equation Q shows that the normal pressure should 
scale with the square of the shear rate. Indeed these scaling laws were measured experimentally in a Couette 
geometry with a fixed volume |I4j . 


2.3. ’’Local” rheology 

Assuming the previous equations govern the flow properties locally, one can predict the velocity profiles 
in other configurations. Let us first consider the inclined plane where the tangential to normal stress 
ratio is fixed by the slope. Using equations to predict the velocity profile, the so-called Bagnold’s 
profile is recovered in quantitative agreement with measurements 0 0 . The second validation was the 
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Figure 2: Schematic of the evolution of the friction coefficient fi (left) and of the volume fraction (p (right) as a function of the 
inertial number 1. 


two-dimensional flow on a heap. The concave velocity profiles are recovered by this rheology if one takes 
into account the influence of the sidewall correctly [TS] (see section . 

To describe 3D flows, such as flows on complex topography or with transverse variations, a tensorial 
form of eq. ([^ was proposed m- The deviatoric part of the stress tensor aij = —PStj + Tiji 

Tij = where I 7 I = (3) 

To write the above equation, the co-linearity of the shear stress and strain rate tensors is assumed 
as suggested in first approximation by theoretical considerations and numerical simulations [ElEQlEI]. A 
closer analysis of the numerical simulations reveals however that there is a non-zero second normal difference 
of roughly 10% of the mean pressure (see [50] and below). Let us note further that equationsandneed 
a closure equation. If the pressure is not constant in time, those equations cannot be used to compute it. A 
simple assumption commonly made in numerical works is to consider the granular flow incompressible. On 
the other hand to capture specific volume-fraction effects, an additional equation is required for its evolution 
(e.g. for immersed dense flows [55]). 

The first experimental test of this 3D law was done with flows on heap between rough walls which 
produce strong transverse shears, concluding that the model captures accurately the transverse free-surface 
velocity profile m- However, the model still predicts a finite thickness instead of the slow, exponentially 
decreasing flow deep in the pile. Despite this limitation, these correct predictions allow to gather the bulk 
flow properties from different configurations in a same theoretical framework. 

2.4- Complex geometries 

The /i(/) rheology was derived from simple flows, and first validated using steady unidirectional flows. 
In front of its successes, researchers have studied more complex dense granular flows: To test further the 
prediction of this rheology, two main approaches were used. A first one is to simulate the granular flow as a 
continuum media using Navier-Stockes equations |231124] . However using this pressure-dependence rheology 
is difficult, mainly because of the divergence of the viscosity close to the arrest. A second approach involves 
using discrete element method (DEM) to access locally to the rheological parameters in the complex flow 
and look for relation between shear stress and strain rate. I present some recent works supporting the ^{I) 
rheology carried on less constrained flows. 

The hrst question that arises is whether this law is still valid for time-dependent flows, i.e. for accelerating 
flows. Studying the starting of unidirectional flows on heap, it was shown that both the time evolution and 
the amplitude can be predicted by this model [25] . A more drastic configuration is the collapse of a granular 
column on a rough plane. This peculiar configuration presents very fast flows that never reach steady states 
with a phase of strong acceleration at the base of the column [55] HT] [55] . Lacaze et al. performed numerical 
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collapses of column by DEM, then analyzed the stress and strain rate tensors. Plotting the local friction 
coefficients as a function of the local inertial number, they obtain a collapse of all their data [2^, although 
they observed a misalignment between tensors of few degrees in the foot (the static part) of the column. 
Recently, using a continuum code, dynamics of the collapse of a 2D column was achieved. The very good 
agreement of shapes and dynamics between the continuum simulations and DEM simulations confirms the 
validity of a frictional rheology. However the relevance of the shape of ^(/) was questioned since a model 
with a constant friction coefficient could produce fairly good results in this case 

Another free surface flow exhibiting spatial inhomogeneity is the rotating drum configuration. Cortet 
et al. analyzed carefully the numerical discrete simulations and found contrasting results: for the fast thin 
flowing layer that accelerates and decelerates exchanging grains with the erodible bottom, the local rheology 
applies |30j . Deep in the drum however, the co-linearity assumption is no longer valid. The researchers 
attributed it to the evolving contacts anisotropy between grains during the rotation or to the transient time 
of arrest. Surprisingly, the invariant form of equation ([^, |t| = /r(/)P, is observed to be valid down to very 
small values of / (10“®). 

Two last works tested the /j,(J) rheology in vertical chute. Eirst, the silo discharge through a small 
aperture was simulated using an incompressible continuum code. Two major features of the silo drainage 
were recovered: the independence on the height of the flow rate and the so-called Beverloo scaling linking the 
flow rate to the outlet diameter |31j . In a second simulation, the drag force acting on a horizontal rod across 
the chute was simulated with a continuum code |32j showing that the scaling with the Eroude number is in 
qualitative agreement with experiments for slow flows. Let us mention that the incompressibility used by 
the numerical solvers is not always a valid assumption, especially near the aperture or downstream the rod 
where the pressure drops to zero. Including the volume fraction variations remains a numerical challenge. 
It is worth noting that the precise shape of the law /r(/) was not thoroughly tested in these continuum 
codes, but only the frictional viscosity (and yield stress) and the saturation of the friction coefficient for 
higher I. Finally, an important remark concerns the nature of the flows mentioned above: The strain rate 
remains moderate in all configurations I < 0.5 (except in some parts of the column collapse), this issue will 
be discussed in the section [5]2] 

2.5. Beyond the average flow 

Up to now, I focused only on the average properties of the granular flows. However, the grain trajectories 
fluctuate around the mean path. The rheology is unable to take them into account, but measurements shows 
evidence of evolution of their amplitude with respect to the shear rate |33j or to the inertial number I [6]. 
The standard deviation of the velocity fluctuations may be a consequence of the local shear rate and simply 
be a passive scalar [5]. Nevertheless, understanding their origin will allow a correct description of the mixing 
by diffusion of the grains during the flow [34) . 

Moreover, correlation lengths on velocities 0 cni or forces (Hill om were measured and were shown 
to increase sharply when approaching the slow dense regime. Having correlations over longer and longer 
distances signs the beginning of observed deviations from the prediction of the rheology. These growing 
length scales are the base of new theoretical developments to better understand the granular properties and 
to extend the existing model toward slow flows. Some of them will be presented in the next section. 

3. ”Non-local” rheology 

Let us first remind that the term non-local rheology is used in this article to enhance the spatial corre¬ 
lations that govern the flow behavior, even if all equations must still be evaluated locally. 

3.1. Evidence of coupling across space 

Most of the granular flows in geophysics or industry are rather thin compared to their lateral extend. It 
is therefore important to compare the prediction of the rheology to the observation. 

The first signs of the failure of the rheology fi{I) can be seen in the velocity profiles in several configu¬ 
rations. The scaling law on inclined plane was shown to deviate from the Bagnold’s scaling as the thickness 
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of the flow decreases. The curvature of the velocity profile decreases and can be inverted [niiiniii]- The 
granular flow exhibits then a transition from the Pouliquen’s flows rules to a non-local behavior [881139j . 

These behaviors become more important close to the jamming transition. Therefore, such effect will be 
more pronounced at the interface between a fast flow and a quasi-static zone. Indeed, free surface flows on 
erodible material (heap or drum) exhibit a creeping region where the velocity profile decays exponentially 
[m m SO] . This feature is neither captured by the model which instead predicts a finite depth [T7] . The 
same results are observed in cylindrical Couette devices where the local friction coefficient is no longer a 
single function of I but depends on the internal diameter in the quasi-static region |41j . Interesting new 
ideas came from split-bottom geometries [421148j , in which the stress map produces wide shear bands in 
slow flows that do not fall into the p-(/) rheology framework |331[^. Derived experiments and simulations 
have evidenced a mechanical coupling across space away from a shear band (e.g. following the movement of 
an intruder far from a shear band) [IHl SH SZl HH] ■ It is worth noting that these previous behaviors become 
visible in the presence of large stress heterogeneities. 

3.2. Modeling flows with high stress heterogeneity 

Correlations have been incorporated in several models to reproduce the non-local effects. One idea is 
to assuming the dissipation involves more grains (arches or columns). For example, assuming the kinetic 
energy is dissipated fully non-locally (over the whole thickness), a linear velocity profile is predicted on 
inclined plane m, Mills et al. integrates viscous stress along force chains [15], Jenkins et al. dissipates 
kinetic energy of column of grains ED). Andreotti extends a single grain model, which already captures the 
essence of the /r(J) rheology, assuming also that collisions occur between columns of grains in the bulk [5T| . 
These situations lead to an increase of momentum exchange and finally modify the velocity profile. While 
these models provide new features, they rest on the existence of the correlations and yet this experimental 
fact is only supported by few numerical evidences [S2| and is still matter of debate [55] . 

In the following paragraphs, I will focus on situations exhibiting both quasi-static and flowing regions. 
When the stress inside the system is highly heterogeneous (silo, cylindrical Couette, shear plane under 
gravity), shear bands appear and conhne the flow. Nevertheless, far from the shear band, the material is 
still agitated by mechanical fluctuations gSJIllllTlIlHI. 

The main advantage of the approaches discussed below is the ability to describe both the quasi-static 
regime and the flowing one, while incorporating the /i(/) rheology. One idea is to assume that the local 
shear rate provides mechanical fluctuations that promote the movement of other grains [54] . The fluctuation 
intensity decreases exponentially with a single characteristic length scale. This model succeeds in predicting 
qualitative behaviors such as shear banding in silo or the change of concavity of the velocity prohle on 
inclined plane. 

Different recent approaches are inspired by works on soft matter |55j . where elasto-plastic events allows 
the grains to flow below the yield stress, and they have been adapted to the rigid nature of grains. At the 
lower order, these models reduce to a second-order spatial derivative equation of an order parameter (e.g. 
fluidity or inertial number) with a diverging length scale approaching the critical yield stress [5611481157] . 
These models based on numerical evidence can predict accurately the velocity profiles in quasi-static regions 
and show indeed a diverging length scale. It is worth noting that the scaling of this length I is close to 
the one derived by analogy to the Prandlt mixing length |2|: Assuming a viscosity scaling with 
equation Q leads to 1{I) « p.{I)/I. 

The remaining question is now to find the physical origin of such correlation lengths. Several candidates 
are identified arches, contact and force network anisotropy, swirling clusters [58] . the coordination number 
[59] . One strategy of ongoing researches to test further the rheology is to study the response to small 
perturbations of an agitated or sheared steady system [sniisii- More details on the rheological properties 
of slow flows will be found in [53] . 

3.3. Flow threshold 

Another signature of the correlation lengths is the existence of flow thresholds [2] , such as the minimum 
height of a flow at a given slope on inclined plane /istop, the minimum flow rate on a pile or the minimum 
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rotating rate in a drum experiment. All these conditions are not captured by the local rheology which 
predicts steady flows at thicknesses as low as possible. 

To address the issue of the arrest of flows, it seems natural to introduce a length scale able to correlate 
the movements of the grains over the whole flow thickness. In the ’’self-activated flow” model [^, the flow 
ceases when the amount of the fluctuation is not enough to sustain it due to the reduce thickness. However 
the hysteretic nature of the arrest and start of the flow is not captured. The starting-angle dependence on the 
thickness of an avalanche has been recently addressed through the analysis of stable granular configurations 
with some success [64] . 

Finally, these thresholds depend not only on the flowing material but also on the roughness of the base 
[TIiSlISSj- Tuning the amount of injected energy at the base in DEM simulations can modify the flow 
threshold and the bottom shear rate, but not the upper part of the velocity profile m- This dependence on 
the nature of the dissipation at the base is difficult to understand if the correlation-length is only governed 
by the bulk behavior. However this length may be defined locally and may depend on the local agitation as 
discussed in the previous section. This issue raises the difficulty of setting correctly the boundary conditions. 
In the next section, the different kinds of interactions between grains and boundaries are presented. 

4. Surface rheology 

Except from some numerical simulations, in every experiments or processes, one has to deal with lateral 
walls, fixed bottom. So besides the bulk rheology, one has to understand the properties induced by these 
different boundaries. The behavior of granular matter contrasts with the one of classical fluid where a zero 
velocity condition is often achieved. Since there is no scale separation between the microscopic length scale 
(the grain) and the typical length of the flow, a change in the roughness of a boundary may dramatically 
modify the flow properties such as the slip velocity [68] . Which property is governed by a given boundary? 
Is it the shear rate, the average slip velocity or the shear stress? 

I review in the following the current knowledge on the influence of boundaries. 

4.I. Lateral confinement 

One of the most underestimated effects in the past years was the effect of the lateral boundaries. Although 
it was noticed early that the presence of smooth sidewalls modifies slightly the flow regime [5S] , understanding 
their effect was the key feature to apply correctly the fi{I) rheology to the heap configuration [TB]. It was 
shown at the first order that a smooth vertical wall adds a frictional force to the momentum balance. Taking 
it into account and assuming a minimum friction coefficient at the liquid/solid interface allows to predict 
first a finite thickness of the confined flow over a static pile and second the increase of the inclination when 
increasing the flow rate as observed experimentally and numerically [ISIIIHIITO]. However the effective 
friction coefficient that one has to use is somehow lower than the microscopic friction coefficient for these 
materials HHum. This means either that the friction is not fully mobilized, for example due to the rotation 
of particles, either that friction forces are not completely aligned with the flow direction. A convincing work 
by Richard et al. showed that below the fast flowing layer on a heap, the effective friction coefficient at the 
wall decreases due to the random motion of the grains in the quasi-solid phase [72] . Moreover they showed 
that the decay of the friction coefficient is linked to the evolution of the density in the very fast flowing 
layer and that the creeping flow corresponds to a glassy state. Thus the effective friction is not an intrinsic 
property of the grains but depends on the microscopic agitation. It is worth noting that the resulting 
characteristic length of the exponential decrease of the velocity also depends on the fast layer properties. It 
is then not obvious if the non-local models described in the previous section can capture this exponential 
tail. 

On the other hand, the use of rough lateral wall to impose a zero-slip velocity seems easier. Jop et al. 
shows that, although this condition is appropriate at the first order to describe confined heap flows, careful 
analysis of the experimental velocity profiles reveal a slight slip at the rough wall m 

The main shear often comes from the bottom and the sidewalls act as a perturbation of the main flow. 
This may explain why simple laws can be used in first approximation here. These puzzling effective friction 
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coefficients or slip velocities have to be understood to address dense granular processes accurately. Although 
boundary conditions can be derived for kinetic theories, I will now focus on the stresses and slip velocities 
from some experimental and numerical works for basal boundaries. 

4-2. Rough bottom surfaces 

First evidences of the influence of the basal roughness on the flow properties are founded in the inclined 
plane configuration. Changing both the size of the glued spheres on the plane and the one of the flowing 
grains modihes the /istop(0) curves [H |73l |65]. The influence of the roughness on the threshold can be 
interpreted by the change of dilatancy [15] . 

Pouliquen first shown that his scaling law still apply: The average velocity scales with the inverse of 
^stop corresponding to a couple grain-bottom. However, simulations of inclined plane have shown that the 
bulk flow properties are not influenced by the boundary conditions m- Moreover contrasting results exist 
in the literature: Shear plane simulations exhibit a rather homogeneous inertial number while inclined plane 
simulations show a strong variation of the inertial number close to the boundary [751166] . It is thus not clear 
if changing the roughness of the bottom only modifies the velocity and shear in a boundary or if it changes 
the whole bulk rheology [55] and so if the link with hstop is a coincidence or not. Non-local rheologies would 
provide some answers. 

The role of the roughness has been thoroughly studied measuring the movement of a single grain down 
on inclined plane. An analytical 2D model captures the same phenomenology than the rheology fj,{I) [51] . 
If disorder is introduced in the position of the glued spheres but keeping the same density, then the same 
mean velocity is recovered I751IIZ1. 

Koval et al. [zaizs] shows a sharp increase of the apparent friction coefficient at the wall when increas¬ 
ing the roughness above 0.4 then a saturation. Moreover, the slip velocity increases when decreasing the 
roughness, in agreement with previous results m- 

Beside the slip velocity, changing the roughness also modifies the structure of the flow. When the ratio 
of the grain diameter at the boundary and in the flow decreases below 0.61, a transition occurs between a 
disordered system to an organized state in compact layers [551 [SI]- Moreover, if the rough bottom is made of 
a regular dense ordered packing, oscillations between layered and disordered systems can be observed [S5] . 

Finally, the microscopic friction coefficient are shown to play a role only for small roughnesses [79l l83] . 
while the increase of the rolling friction coefficient can trigger a transition between a sliding mode to a rolling 
one [SI]. 

4-3. Smooth bottom surfaces 

It is expected that changing a rough bottom plane to a smooth one induces a slip velocity at the 
interface. However, the bulk properties such as the shear rate and volume fraction do not seems to be 
affected in simulations [S3 IS]- Let us then focus on the interface properties. 

For very low values of wall-particle friction coefficient g,ujp, the material slips as a solid block. Increasing 
above a critical value leads to the development of a shear layer [S3] . A hrst surprising result is that this 
critical value is lower than the particle-particle friction coefficient. The second results is that the effective 
friction coefficient (ratio of tangential stress to normal stress) drops significantly mimicing the static/dynamic 
transition for classical solid friction [86] . In a smooth shear plane conhguration, they identified a boundary 
inertial number which governs the friction coefficient and three different regimes are observed: a linear 
velocity profile above an inertial number of 0.01, two symmetrical sheared layers at intermediate velocities 
and a single shear band regime for very low velocities. 

Even though the link with the bulk inertial number is still lacking, an interesting step has been made to 
relate the slip velocity to the physical properties. Artoni et at studied the movement of a grain submitted to 
a stochastic force from an imaginary bulk and inferred a scaling law for the slip velocity [87] . This empirical 
scaling has been confirmed analyzing flows of angular grains down frictional smooth plane numerically [85] . 
The slip velocity is governed by the ratio of the effective wall friction over the wall-particle friction coefficient 
Heff /IJ'wp- It is worth noting that their law can be interpreted as a hydrodynamic slip length that diverges 
when going closer to 1 in the similar way that the diverging correlation lengths presented in the non-local 
continuum models. 


Let us finally notice that a smooth frictional wall cannot model a rough plane: The two effective friction 
coefficients coincide only in the limiting cases of very large or very small roughness |86j . This result reminds 
of the complex role of the roughness already mentioned in the flow threshold section. 

4.4- Erodible interface 

In the case of a flow over an erodible bottom, different approaches are possible. As shown previously 
[nins], if the ground is composed of the same material then the bulk rheology could be used to find a first 
estimate of the evolution of the interface with time. Two difficulties arise. Since these flows are unsteady, 
there is no simple law linking the flow properties and the position of the interface and one has to compute 
it numerically. Moreover this ’’interface” separates fast and quasi-static flows, then the non-local effects of 
slow flows are more pronounced (see section]^. 

In the following I briefly discuss an alternative approach to overcome this issue integrating the momentum 
and mass equations over the height. When considering avalanches, an additional equation is required to 
couple the interface and flow properties evolutions. A first approach is to set the amplitude of exchange rate 
with the distance to a neutral inclination [88l [89l [90] . Another idea is to impose a property of the flow (the 
velocity profile or the basal shear stress) then moving the interface modifies both the height of the flow and 
its average velocity Eiiiniisa. Such method is useful for example to describe the wave instability observed 
at high flow rate |94| or to study the mobilization of a resting granular bed during the transient regime of 
an avalanche |95| . However the quantitative link with the bulk rheology is still lacking and required more 
effort on deriving analytical models. The case of avalanches on very shallow flows [Miinzj can be described 
using a solid/liquid mixture model whose ratio is governed by an a Landau equation on the volume fraction 
[98) . It could be very fruitful to try to model them using the recent theoretical developments. 

A last case of importance for the geophysical applications is when the ground has a given cohesion. Not 
only the average rheology at the interface is crucial to predict the velocity profile, the ground behaving as a 
rough solid bottom, but also the fluctuations of forces are important to understand the rate of the erosion 
process [99]. Understanding the physical background of these fluctuations, the impacts of grains, will also 
provide answers on the origin of the rheological law [10011101) . 

5. Discussion and conclusions 

5.1. Slow deformations 

So far we have seen that the rheology /i(/) presents a minimal relation to describe a large variety of dry 
granular flows. Moreover attempts have been made to extend the domain of application or to include more 
mechanical behaviors: e.g. to combine this rheology to the hstop |54) . to identify a relevant parameter for a 
non-local rheology able to explain creeping flows far from high sheared areas iMiiiainz]- 

The remaining open question is the physical origin of these non-local behaviors. Obviously the mi¬ 
crostructure is important. The contact network and force anisotropies are shown to reflect the evolution of 
the friction coefficient even at high inertial number I in shear plane flows [102) . However, it is not clear 
whether such consideration could help explaining non-local effects: Which parameter governs the dynamics 
in this case? Minute variations of the volume fraction, assumed to be a passive parameter linked to I, could 
play a role in slow flows. The importance of the coordination number in the finite size effects on inclined 
plane (hstart) has been also pointed out [64]. 

Although such approaches present a substantial improvement, some properties of the slow granular flows 
would still not be captured. In the inertial regime only slight normal stress differences have been identified 
m- However, evidences of these effects can be found in slow flow experiments, such as the levee formation 
[inaEHiisi] or free-surface transverse curvature uns]. These effects can be casted into the local 3D rheology 
[105) and the depth-average predictions are in good agreement with the experimental surface proflles. Further 
tests of their model are now required on 3D flows. 

Finally, let us mention theoretical developments of more general mechanical models for plastic deforma¬ 
tions in granular materials [106) that could serve as guide to identify the key parameters. 
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5.2. Fast flows 

The other limit of the rheology is the fast flow regime when part of the system becomes dilute. 
Transition toward salting grains are observed when the shear to normal stress ratio exceeds the limit fj.g 
[ziiniin]. Numerical results show that beyond a critical value, an unstable branch of the friction coefficient 
exists: The friction coefficient starts decreasing with I for greater values than 0.8, corresponding to the 
beginning of the dilute regime fTU7llilig|. The precise value of the transition depends on the softness of 
particles [108j . These finding are in agreement with the low friction coefficient found in the single-grain 
model when the grain accelerates on a bumpy surface |51j and with the intermediate friction coefficient 
for numerical supported flows ung. The description of very dilute flows relies on the physics of binary 
collisions and predicts indeed a decrease of the friction coefficient with the inertial number. However, dense 
stationary flows are still observed experimentally at high inclination, either in a confined configuration [13] 
or for shallow flows on plane m- These flows strongly challenge the rheology: In between the dense 
regime at moderate inclinations and the complete dilute regime, the flow remains relatively dense and large 
deviations from the predictions of the /i(/) rheology are measured with higher friction coefficients [llOlllll] . 
Aside from the stationary flows, instabilities often arise from the spatial variation of density and these new 
patterns, out of reach of the fi{I) rheology, strongly narrow the range of application of the rheology and can 
serve to test the kinetic models [miniiiiniiisi. One crucial issue would be to identify the nature of the 
dissipations of energy in such flows: e.g. inelastic collisions, secondary Hows, ’’turbulent” flows. 

To adapt the kinetic theories to the dense regime, researchers have relaxed the assumptions that the 
flow is dominated by binary collisions, since neglecting long lasting frictional contact may introduce large 
errors m- Extensions to adapt the kinetic models to the dense regime provide finally a correct diverging 
thickness of hgtop for low inclination |115l I116j and different features of inclined plane and heap flows can be 
well modeled [nil EH]- However the flow is still described as the spatial juxtaposition of different regimes 
(base, dense core, diffuse free surface). 

More details on the recent developments of rheological models based on kinetic theories will be found in 

m- 

5.3. Complex interactions 

The main focus of this article was on the rheology of monodisperse dry spherical grains. What do these 
behaviors become if the interaction law between grains changes? As soon as we introduced a new physical 
ingredient, / is no more the single dimensionless parameter in the system. Several works have tried to 
understand the rheology of polysdisperse systems, cohesive grains or non-spherical ones. 

Pouring grains of different density or size on an inclined plane leads generally to a segregation phe¬ 
nomenon. Even though the mechanisms of the segregation are still matter of debate |120j , recent studies on 
binary dense flows show the validity of the local rheology in steady states. The friction coefficient is shown 
to follow the same dependence with the inertial number provided that the definition of I is modified. One 
should replace the grains diameter by the average one |12111122j and the particle density by its mean value 

[HH]. 

Also the nature of the interaction can be tuned. Attractive forces can arise when decreasing size or adding 
capillary bridges. The rheology of cohesive grains follows the same trend of increasing friction coefficient 
with I, however the cohesion increases globally its values |124LI125] . This effect is linked to the modification 
of the microstructure with an unexpected role of the distant attractive forces |126j . 

When the grains are immersed in a viscous fluid, the same frictional behavior is observed and can be 
described provided to change the microscopic inertial time scale by a viscous one ina- Aside the dynamical 
effect of volume fraction evolution, the very same theoretical framework can be applied then to dry grains 
and dense suspensions |1281112911130j . 

A final parameter is the shape of the particle. Irregular grains, like sand, deviate from the prediction 
of the simple rheology. Elongate grains are shown to orientate with the shear direction |131j . Although 
influence on packing fraction and on slow shear regime of non-spherical particles has been investigated, the 
dense flow of such particles deserves more attention in the future. 
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B.Jj.. Conclusion 

I restricted the scope of this article to dense flows. I presented their mechanical properties in the bulk as 
well as near boundaries together with the most recent theories. The average features of the dense granular 
flows are now well understood. Several models have been built progressively to explain them. The rather 
useful model /i(/), presented in the GdR MiDi’s article ten years ago, can predict some of the basic behaviors 
of the dense granular flow at moderate strain rate. Since then the number of configurations of greater degrees 
of complexity that take place within the same theoretical framework has continued to increase. However, 
some common geophysical flows or some industrial applications, such as fast dense avalanches or high-shear 
mixing processes, cannot be modeled accurately by this rheology. This is part of its current limits that I 
have shown and that still require other theories to be understood. 

To conclude, let us mention other open questions. When deforming a granular media very slowly, for 
example in a biaxial test, the material does not flow homogeneously but present instead localization of the 
shear rate. Those shear bands qualitatively differ from the one developing close to a wall or in the split 
bottom geometry. Up to now, they are not described by the same model as dense flows. The question is then 
to try to link the behavior at high and low strain rates. How can the most recent theoretical developments 
describe time-dependent flows where the anisotropy of contact network may be enhanced, such as oscillatory 
flows or the time-evolution of forces around intruder passing through his own wake |132) ? Finally, what are 
the physical origins of the efficient models presented is a crucial question. Underlying the non-local effects, 
I hope these theories will motivate interesting works to enlighten the relevant parameters. 
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